Entropy-driven pumping in zeolites and ion channels 
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When two binary solutions are separated by a permeable barrier, the individual species typically 
diffuse and mix, dissipating their chemical potential gradients. However, we use model lattice sim- 
ulations to show that single-file molecular-sized channels (such biomembrane channels and zeolites) 
can exhibit diffusional pumping, where one type of particle uses its entropy of mixing to drive an- 
other up its chemical potential gradient. Quantitative analyses of rates and efficiencies of transport 
are plotted as functions of transmembrane potential, pore length, and particle-pore interactions. 
Our results qualitatively explain recent measurements of "negative" osmosis and suggest new, more 
systematic experiments, particularly in zeolite transport systems. 
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Particle transport across microscopic pores is a crucial 
intermediate step in almost all biological and chemical 
engineering processes. Separation, catalysis, and drug re- 
lease all rely on controlled transport through microscopic 
channels such as zeolites |^ . Biological examples include 
integral membrane proteins which traverse hydrophobic, 
highly impermeable lipid membranes and participate in 
small molecule transport Q . 

Water transport is usually assumed to be driven by 
osmotic pressure differences across pores Solute per- 
meability may contribute a counterflowing solute current, 
which often reduces the magnitude of total volume flux 
iQ,^. However, we show that for large enough asym- 
metry between solvent-pore and solute-pore interactions, 
volume transfer can occur in a direction opposite to that 
expected from simple osmosis Here, one species of 
the otherwise counterflowing pair wins and forces the 
second to be coflowing. Analysis of a physically moti- 
vated kinetic model of single-file pores shows that they 
can also function as symport pumps [Q under certain 
physiological conditions, where for example B(=solutes, 
ions) pumps A(=water) up its electrochemical potential 
gradient. Recent osmosis experiments in various sys- 
tems have indeed revealed negative reflection co- 
efficients, i.e., volume flow opposite that expected from 
standard osmosis. 

Since transmembrane pores have typical radii ^ 1 — SA, 
we use a one-dimensional exclusion model to de- 

scribe single-file particle flows. A single-file chain is di- 
vided into N sections labeled i, each of length £ > molec- 
ular diameters of A and B, and containing at most one 
particle of either type. The A(B) occupation at site i 
is defined by Ti,T^ £ {0,1}. All parameters associated 
with A(B) will carry unprimed(primed) notation. The 
probability per unit time that an A(B) particle enters 
from the left i/ site i = 1 is unoccupied is axL{ct'x'L)i 
where Xj x' are the reservoir mole fractions of the partic- 
ular types of particles that enter. Similarly, the entrance 
rates at i = iV are denoted by SxRiS'xR)- We assume 



= 1 - XL > Xfl = 1 - XR^ as depicted in Fig 1(a). The 
A(B) exit rates to the right(left) are denoted by /3, 7 and 
/?', 7'. In the chain interior. A- type particles move to the 
right(left) with probability per unit time p{q) only if the 
adjacent site is unoccupied. B-particles hop left(right) 
under the same conditions with probability p'{q'). 
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FIG. 1. (a) A two-species ID pores representing e.g. zeo- 
lites or biological transport channels is divided into sections 
of length I. The kinetic rate constants for A (unprimed) and 
B (primed) particles are indicated, (b) and (c) depict energy 
barriers for A (open circles) and B (filled circles) particles. 
Assume A is uncharged (AV = Vr — Vl = Q), while B may 
be singly charged and acted upon by a transmembrane po- 
tential A1/' = Vr-V[ <G. 

Rigorous results exist for the two species asymmetric 
exclusion process under specific conditions | [lO| |. How- 
ever, currents for N < A can be easily found analytically 
for general parameter sets by solving linear rate equa- 
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tions coupling all possible particle occupancy configura- 
tions H^]. For > 3, lattice simulations describing the 
pore dynamics are implemented by choosing a site i and 
finding the instantaneous current between i and i + 1: 
J, = pT,{l - T,+i)(l - T^+i) - gTi+i{l - r,)(l - t[) e 
{0, ±1}, where p ~ 1 with probability pdt, and zero with 
probability (1 — p)dt as dt 0. Analogous expressions 
hold for the distribution of q. The occupations t^jt/ are 
then updated and the next particle randomly chosen. 
An analogous rule holds for J' {(,')■ Boundary site ki- 
netics are correspondingly determined by a, S, a', (5', e.g. 
Jn = /3rjv — (5(1 — TAr)(l — T^). Typically 10^"-^-^ steps are 
needed for J, J' to converge to their steady state values, J 
and J'. We verified all results by comparing simulations 
for both a' = 5' ^ l^V = AF' = with exact resuhs Q 
and general kinetic parameters with numerical results for 
iV < 3 

Since particle fluxes across microscopic channels are 
< lO^/s [||, typical fluid particles suffer > 0{W) 
"collisions" and relax to local thermodynamic equilib- 
rium (LTE) while traversing the pore. Kinetic rates, 

} = {a, /3, 7, (5, p, under LTE take Arrhenius forms: 



p~q^ {vt/C) exp{~Ep/kBT) 
/3 ~ 7 « {vt/£) eM-Ep/kBT) 
a, i5 ~ ao exjp{-Ea/kBT), 



(1) 



with 

analogous expressions for {£_'} = {a', /?', 7', S',p', q'}. We 
have for simplicity assumed microscopically symmetric 
pores and equal hydrostatic pressures {i.e. (3 = = 5 
and P' — 7', a' — 5' when transmembrane potentials 
/S.V = and ISV' — 0, respectively) and also that within 
a narrow pore, particles equilibrate with the pore inte- 
rior and relaxes momenta much faster than particle po- 
sitions, implying, in the absence of external potentials, 
j3 ~ g, y ~ q' . The hopping rates f3,j,p given in (|^) 
represent ballistic travel times (thermal velocity vt di- 
vided by section length), weighted by internal interac- 
tion energies Ep. The entrance energy dependent factors 
(a, 6) and (a', <?'), when multiplied by the relevant num- 
ber fraction of entering particles, {xl,Xr) and {x'l'X'r), 
respectively, define entrance rates into empty boundary 
sites. The precise values of the prefactors ao, will also 
depend on equilibrium parameters in the reservoirs such 
as temperature, total number density, and the effective 
area of the pore mouths. The energy barriers experi- 
enced during single particle transport are enumerated in 
|l|(b) and (c) and may include external potentials. If say, 
the B-particles have charge ze and are acted upon by 
a pondermotive transmembrane potential AV ^ 0, the 
energy barriers are shifted: E'^ ^ E'^ ^ + kgTv' /2 and 
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kBTv'/2, where v' = AzeV /{N + l)kBT. 



Thus, the kinetic rates when v' ^ Q are 



(a',/3',p') 
(7',^',9')- 



E(a',/3',p')e-"'/2 
(7',<5',p')e+"'/' 



(2) 



When A(B) is uncharged(charged), J, J' are computed 
using the parameters = (a'^, P'v,7v,K,PviQ'v)^ 

and x'l,r- 

Under isobaric, isothermal conditions, AE{Pl = 
Pr) = Ea + Ef3 — Ej — Es = 0, and since pressure fluc- 
tuations in liquid mixtures equilibrate much faster than 
concentration fluctuations, the total enthalpy change per 
A(B) particle translocated is AH{AH') ~ AV{AV'). 
The efficiency of using species j to pump k can be de- 
fined as the ratio of the average free energy gained by k 
to the free energy lost by j: 



r^,k = [l-e{.PA^i^)\ 0(J^AAi^) 



(3) 



The Heaviside functions represent the definition that 
fiow is considered useful work only when j and k are 
coflowing. Using the entropy of mixing per particle 
AS" = fcBln(x^/x^), the Gibbs free energy change per 
particle, A^' = AH' - TAS', is 



A^' ~ fesTln 



x'l 



AV, 



(4) 



with an analogous expression for Afi associated with the 
transport of an A-type particle. When concentration 
changes, are not too large, higher interaction terms con- 
tributing to AH, AH' (and hence A/i, A/i') can be ne- 
glected. These correction terms (which are higher order 
polynomials in Ax, A^') can be straightforwardly incor- 
porated by independently measuring bulk liquid heats 
of mixing. For concreteness, we assume species j = B 
(charged) is used to pump k = A (uncharged) and that 
for liquids under ambient conditions vx/i ^ lps~^. Upon 



setting the time scale dt 
1. 
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FIG. 2. 10r?(iV = 3) (exact, solid curves), 10r?(A^ = 5), 
and 5ri{N — 10) as functions of AV' . The cor- 
responding parameters (for a B-attracting pore) are 
l3 = f = p = p' = q = q' = 0.01 (dt = lOfs), and 
P' — "f' = 0.003. The intrinsic entrance rates here are 
Of = 5 = q' = 5' = 0.01. The smaller f3' < (3,p,p' cor- 
respond to a pore that is more attractive to B than to A. 
Here, x'l/x'r = (1 - Xl)/{1 ~ Xr) = 0.0036/0.0018, corre- 
spond to a 200mM/100mM aqueous B solution in {L)/{R). 
The computed currents and efficiencies are low and (numeri- 
cally) noisy due to low (physiological) concentrations x\ Bt,nd 
entrance probabilities a'x' ■ 
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Fig. H shows currents and efficiency as functions of a 
transmembrane potential difference Al/' ^ 0, for var- 
ious length pores. Note that under the physiological 
conditions considered, the efficiencies are small (< 1%) 
for pores of molecular lengths. A small |A/i'| would in- 
crease efficiency via transpore energetics; however, for 
too large a AV', the B-particles are driven back against 
their number gradient, and useful work precipitously van- 
ishes. This occurs most easily for short channels where 
internal A-B interactions are rare; here AV' < is re- 
quired for B to drive A uphill. The efhciency is non- 
monotonic and has a maximum as AV' is varied for fixed 

U,nandxi.,fl. 

For AV' — > — oo, an asymptotic form for the efficiency 
can be found by assuming the B particles never hop 
against the potential gradient. 
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Although J', which is dissipating down its electro- 
chemical potential, generally decreases for longer pores, 
for smaU AV, J{N = 10) > J{N = 5) > J{N = 3) be- 
cause pumping is more efficient since there is less likeli- 
hood that a B particle can drift through without pushing 
out all the A particles ahead of it. 
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FIG. 3. Dependence of (a), currents and (b). efficiency 
on pore length N = L/l for zeAV = 0.0, O.SfcsT. The pa- 
rameters {£,, ^'} used are the same as those in Fig. except 
x'l/x'r = 0.18/0.09 corresponding to a 10M/5M aqueous so- 
lution. Although J, J' generally decrease with A'', efficiency 
increases to an asymptotic value. 

The particle fluctuations across the pore that allow A 
and B to simply dissipate their chemical potentials be- 
come rarer as membrane thickness or N increases, en- 
hancing efficiency as shown in Fig. |3[ Henceforth, unless 
otherwise indicated, we treat only entropic driving, i.e., 
AV = AV = 0. Note that r]{N = Q, 1) = is exact for 
all parameters since = corresponds to an infinitely 



thin, noninteracting membrane, and an analytic solution 
for J{N = 1) oc aPixR - x'l) 0- When A/z > 0, 
x'b. x'l J < 0, regardless of whether or not solute 
enters or passes through the single-site pore. This is ex- 
pected since A and B never interact within the single-site 
pore for B to be able to "ratchet" A through. For larger 
N however, we find an asymptotic maximal efficiency 
defined by the kinetic parameters {C)^'}- As A^ — > oo, 
the large fluctuations required for net particle transport 
will transfer a constant ratio of A and B particles. This 
asymptotic efficiency can be tuned by judiciously select- 
ing {Cj'C'} which yield the desired N ^ oo performance. 
Note that for AV = 0.3, the 77(A^ oo) limit is larger, 
but for short pores requires larger A^ both J, J' > and 
pumping to take effect, consistent with the results in Fig. 
2. 
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FIG. 4. (a). Currents and (b). efficiencies as functions 
of B-pore binding (~ 1//3')- The solid curves correspond 
to a slightly A-attracting pore where a — 0.01, /3 = 0.005. 
The dashed curves correspond to an A-repelling pore where 
a = 0.005, /3 = 0.01. 

Finally, performance can also be controlled by micro- 
scopic pore-molecule interactions. For example, recent 
simulations have demonstrated dynamic separation fac- 
tors between species with different pore binding char- 
acteristics We suggest here that similar molecular 
considerations can be used to design pores ||l5|| that op- 
erate as thermodynamic pumps, without moving parts, 
complex cooperative binding mechanisms [0,^, or ex- 
ternal driving forces. Fig. ^ shows the effect of varying 
B-pore binding ^ 1 //?' for two different values of A-pore 
affinity a/ (3. 

Pores that do not strongly attract (larger ||l^) , 
or that repel solutes (smaller a') have lower B particle 
average occupations r', and allow solvent (A-particles) 
to more likely pass from (i?) to (L), decreasing efficiency. 
However, at very small P' (highly attracting pores), r' in- 
creases to the point where it fouls the pore and decreases 
J relative to J', also decreasing efficiency. Thus, there ex- 
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ists an intermediate B-pore affinity, which yields a maxi- 
mum efficiency. Depending on A-pore affinity (dashed vs. 
solid curves), currents and efficiencies can be simultane- 
ously increased {e.g. by making pores slightly A-repelling 
when 0.0028 < P' < 0.004). These behavior are a conse- 
quence of the intrinsic nonlinearities and are not present 
in the phenomenological linear Onsager limit 

The above results are consistent with, and provide a 
microscopic theoretical framework for recent anomalous 
or "negative osmosis" experiments ^. These measure- 
ments show negative osmosis only when the chemical 
structure of cation exchange membranes was modified 
by adding — CH2— functional groups. Thus presumably 
makes the pores smaller and enhances the likelihood of 
finite particle size exclusion. More definitively, negative 
osmosis was measured for a window of fixed membrane 
charge (— CH2— SO3 ^) density. Tuning these fixed coun- 
terion charges is equivalent to tuning (3' for positively 
charged solutes (solutes which exhibited anomalous os- 
mosis were Ca+^ Ba+2, Sr+2, and not H+, Na+) §. 
Figure ^ shows a maximal currents and pumping effi- 
ciency as a function of f3' and matches these experimental 
findings. 

There is also evidence for biological manifestations of 
diffusional pumping Recent experiments show 

that water transport is coupled to Na^-glucose [0 and 
KCl Q transport. The negative osmotic reflection coef- 
ficient measurements across Necturus gallbladder epithc- 
lia [Q in particular have eluded explanation, although it 
is conjectured that separate KCl and water transporters 
must be near each other in the membrane and coupled 
Q. However, we have demonstrated how a single sim- 
ple pore can exchange free energy (even in the absence 
of direct forces such as transmembrane potentials) be- 
tween two components and utilize entropy to perform 
work. This implies that under certain conditions, com- 
mon membrane channels can mimic symport pumps 
which are conventionally thought of as more complicated 
shutthng proteins |p|,[7|,^ . 

Biological cell membrane channels have sizes that limit 
diffusional pumping efficiencies, particularly at physio- 
logical solute concentrations (Fig. 2). However, the fiow 
rates achievable by simple pores are higher than those of 
shuttle enzymes and may be a viable mechanism in cel- 
lular volume control, or whenever high fiuxes are desired. 
The ubiquity of membrane channels that conduct water 
1^,^ , and are leaky to certain solutes , suggests the 
mechanisms presented should be considered when inter- 
preting "negative" osmosis and coupled transport exper- 
iments. 

We have not treated "slippage," or incomplete coupling 
processes defined by o» ^ #0 which may occur in wider 
channels and decrease efficiencies. Similarly, attractive 
interactions between A- A, A-B, and B-B can also affect 
performance either way. More accurate molecular dy- 



namics or Monte-Carlo simulations may reveal further 
details of diffusion pumping. Systematic measurements, 
especially on more controllable artificial membrane sys- 
tems where a wider range of the parameters we have con- 
sidered can be explored, may eventually reveal secrets of 
more complicated bioenergetics. 
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